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Abstract 

Several boson exchange potentials, describing the NN interaction Ti < 300 MeV with high 
quality, are extended in their range of applicability as NN optical models with complex local or 
separable potentials in r-space or as complex boundary condition models. We determine in this 
work the separable potential strengths or boundary conditions on the background of the Paris, 
Nijmegen-I, Nijmegen-II, Reid93, AV18 and inversion potentials. Other hadronic systems, 7rN, KN 
and 7T7T, are studied with the same token. We use the latest phase shift analyzes SP00, SM00 and 
FA00 by Arndt et al. as input and thus extent the mentioned potential models from 300 MeV 
to 3 GeV . The imaginary parts of the interaction account for loss of flux into direct or resonant 
production processes. For a study of resonances and absorption the partial waves wave functions 
with physical boundary conditions are calculated. We display the energy and radial dependences 
of flux losses and radial probabilities. The results lend quantitative support for the established 
mental image of intermediate elementary particle formation in the spirit of fusion. 



* To appear in the Procedings of the International Workshop on Resonances in Few-Body Systems, 4-8 

September 2000, Sarospatak Hungary, Few-Body Systems, Supplement, Springer Verlag 
^Electronic address: geramb@uni-hamburg.de 



1 



I. INTRODUCTION 



From the recently fitted high quality boson exchange models their is none which fits 
the latest NN phase shift analyzes, which today extents to 3 GeV[0]. Several theoretical 
attempts, which explicitly included the Delta and other nucleonic resonances, remained qual- 
itative and do not meet the required high standard of reproduction of phase shifts or on-shell 
t-matrices, as is needed in their use for nucleon-nucleus reaction studies. Potentials are still 
very convenient to extent the two nucleon on-shell t-matrix into the off-shell domain which 
enters in few and many body reaction calculations. It is now well established that all the 
high quality NN potentials have a very similar off-shell continuation which is practically the 
same as implied by local inversion potentials in r-space. Repeatedly occur claims, regarding 
the NN t-matrix off-shell, of sensitivity in few body reactions. Practically, these claims end 
mostly in clear failures of precise on-shell reproductions of t-matrices. We followed such 
studies in pp bremstrahlung and microscopic optical model calculations at low and medium 
energy and found unanimous support for need of a precise on-shell reproduction with little 
effects from far off-shell t-matrices. This is what our approach satisfies. 

Furthermore, the 7rN boson exchange models and older phenomenological separable po- 
tentials describe the Delta resonance rather as a molecular system, composed of a well 
separated tt and nucleon, and not as a compactly fused elementary particle. Contrary to 
these findings support elastic scattering inversion potentials in the [3/2,3/2] channel the 
intermediate formation of a fused 7rN system, the Delta resonance. 

Since, already at low energy, production and in-elasticity is important for many elemen- 
tary systems we develop optical potentials for NN, 7rN, tttc and KN scattering. The details 
of how these potentials are calculated with a stringent link to experimental phase shifts is 
contained in Sect. 2. The complex optical model potentials are then used to calculate wave 
functions, radial probabilities and the radial losses of flux from the continuity equation. 

What we expect from these studies is the support of existing mental images of formation 
and fusion of elementary particles. The transition into a combined object, viewed in terms 
of the relative distance between the centers of particles, is to be signaled by a concentration 
of loss of flux from the entrance channel and thus disappearance into the QCD sector and 
deeper down into complex particle production processes. We expect a typical distance of 
0.3 — 0.6 fm with the implication that the charge form factors overlap significantly, 80-95%. 
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Ad hoc, this picture is not supported for 7rN boson exchange potentials, at least in studies 
of the A resonance. To illustrate these findings we show samples in Sect. 3. 

The elastic partial wave NN phase shifts are now available from Arndt et al. for energies < 
3 GeV. Caught from the spirit of inverse scattering, we generate a quantitative optical model 
for any of these NN partial waves with a versatile program which allows studies of complex 
separable potentials on the background of the existing NN potentials like Paris, Nijmegen 
and our own inversions. The described procedure, Sect 2, is the inverse scattering branch 
of soft obstacle identification with elastic scattering data and their partial wave phase shift 
decomposition respectively. With the experimentally determined phase shifts and thus the 
on-shell t-matrix, an inversion algorithm is a mean to continue these on-shell values into the 
off-shell domain within classes of potentials. The class is limited in the present calculations to 
be a range of well established, real valued in r-space, OBEP (Paris, Nijmegen-II, Nijmegen- 
I, Reid-93, AV18, ESC-96) or inversion reference potentials and separable potentials with 
complex strengths. The separable potential strengths are determined with an inversion 
algorithm which guarantees full agreement with any experimental phase shift in any partial 
wave and energy. A versatile range of options for the separable potential form factors permit 
also studies of boundary condition models which suit distinctions between hadronic and QCD 
domains. 



II. THE OPTICAL MODEL POTENTIALS 



NN scattering is formally described by the Bethe-Salpeter equation 



M = V + VGM . 



(1) 



It serves as ansatz for the Blankenbecler-Sugar reduction which is obtained from the integral 
equation (JTJ) in terms of four-momenta 



M(q', q; P) = V(q', q; P) + J d 4 k V(q', k; P) Q{k' P) M(k, q- P) 



(2) 



using the propagator 



G(k;P) 



(27T) 



[\P + k) 2 — M 2 + is 



(!) r 



\ P+ tf+M 
\P + k) 2 - M 2 +ie 



(2) 



(3) 



3 



The superscripts refer to the nucleoli (1) and (2) respectively, and in the CM system, P = 
(y/s, 0), with total energy y/s. The reduction of the propagator Q uses the covariant form 

where only the positive energy projector is used. Finally, a three-dimensional equation is 
obtained from 

AW.qHVW.O + Z gWklf ^V q) ■ (5) 

Taking matrix elements with only positive energy spinors yields the minimum relativity form 

/d z k M 2 1 

W V(q ' k) ^?^T¥ T(kiq) - (6) 

With the substitutions 

n<™-(£)W<)\ v«>Hih<'<t (7) 

we obtain a Lippmann-Schwinger equation for the t-matrix 

/H 3 k M 
^(q / ,k) q2 _ k2 + ^ T(k,q). (8) 

The relation between t-matrix and potential, or free wave function and scattering wave 
function, T$ = V^f gives the Schrodinger wave equation 

[-A + MV(r)-k 2 ]ilj(r,k)=0 , (9) 

where the reduced mass 

nr nr m\ + wi2 

must be used for consistency with reference potentials. A careful and consistent treatment 
of (M/E) factors in the transformation are important for use of t-matrices in few and many 
body calculations where relativity is of importance. Minimal relativity enters only in the 
calculation of k 2 , 

M 2 2 = (mi + m 2 ) 2 + 2m 2 T Lab = (^J k 2 + m 2 + xj k 2 + m%J , (11) 

while the relative momentum in CM system is 

^ = ml{Tl ab + 2m l T Lab ) ^ 
(mi + m 2 ) 2 + 2m 2 T Lab 

which, for equal masses, reduces to k 2 = s/A — m 2 . 



A. Algorithm for Separable Optical and Boundary Condition Models 



We distinguish three Hamiltonians: i.) A reference Hamiltonian Hq := T + Vq with 



S-matrix in the elastic channel), of which any quantity, such as wave functions and phase 
shifts, is readily computed; ii.) A projected Hamiltonian Hpp := PHP, with P + Q = 1 and 
=1 >< (fiil- The Q-space functions \<j>i{r) > are doorway states describing the 
QCD entrance sector from finite size nucleons, mesons and possible other particles. Meson 
creation and annihilation occur only in Q-space and deeper down, iii.) Full optical model 
Hamiltonian TC := T+Vq + V(t, r'\ Isj, E), with scattering solutions := \l/ + (r, k, E). They 
shall match asymptotically continuous energy fit solution of elastic channel S-matrices. In 
the present case, the experimental complex phase shifts are taken from SAIDp|. The refer- 
ence potential Vq and V(r, r'; Isj, E) are to be specified in detail with any of the applications. 

B. Towards Full Potential Model 

Using the scattering equations, (E — H )tp = and (E — H PP )x(jp = 



elastic scattering solutions ipQ := ipo(r, k, E), with real phase shifts (equivalent to a unitary 



(E — Hpp — Hq P — H P Q — H Q Q)ijjp = —HQpijjp 



(13) 



or as Lippmann-Schwinger equation, 



1 



fpp 



(E+ - H ) 



Hqp^p 




(14) 
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Orthogonality PQ = QP = yields 



< i\ip P >= =< i\ip > - < i\G + H QP \ip P >, 



(15) 



and 



N 



< ]\H QP \^p >= J2i< >} 1 < #o > • 



(16) 



The P-space solutions, in terms of ipo, are 




(17) 
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with a separable potential 



> <j\G+\i > < Jl = l? > Xlj < Jl =: A ^ r,r ') (18) 



and 

AT 



^p = ^o-^G + A^ (19) 



The Q-space definition gives a unique value for the separable potential Ay(r, r') and its 
strengths \j{lsj,E). The result looks like a first order Born approximation but is exact. 
However, in the next step we allow the strength matrix 

Ay = {< n\G + \m >}" 1 (20) 

to vary with the implication to reproduce asymptotically the experimental phase shifts of 
the full Hamiltonian 

Hpp^H. (21) 
In case of non unitary S-matrices a complex optical model is the result of matching 

~ ~ ^(-h- + h + S(k)). (22) 



From 



I/^i?)^) - S (k)) = ^ > \ tJ < ^ > (23) 



we obtain the strengths Ay. Next, the separable potential is transformed into 

V ^- X (T^FA) < 24) 

for which the Lippmann-Schwinger equation 

= + G+V# w . (25) 

holds. 



C. Technical Details 

The radial wave functions (we identify the complete channel set of quantum numbers by 
a) of the reference potential 

<M)=(^^ + fK(r,«) 



V{{r,a) V k 2 



1 + 2V b (r, a) J 1 + 2V b (r, a) 



u a (r,k) (26) 



are numerically integrated for single and coupled channels by Numerov. The potentials 
V a , Vb, V b ' are available for the quantum numbers (LSJT) for the Paris, Nijmegen, Argonne 
and our own inversion r-space potentials Paris, Nijl, ESC96 all V b ^ 0, and Nij2, Reid93, 
AV18, HHinv all V& = 0. The physical solutions are matched asymptotically to Riccati 
Hankel functions 

<(r,k)~ l -(h-(r,k)-h+(r,k)S° a (k)) (27) 



and normalized 



^+(r,fc) = ^ P (r,fc) . (28) 
y/l + 2V b (r,a) 



The outgoing wave Jost solution 

j:(r,k)~h+(r : k) (29) 

is calculated by the same token as the physical solution and they are used in the reference 
potential Green's function 

\-m\i> + Ar,k)J+{r>,k), r<r> 
G+(r,r',k) = { _ (30) 

l-(f) lJ+(r,kW+(r>,k), r>r> 

Therewith, at the asymptotic matching radius r = R, the strengths matrix \ a (k) is calcu- 
lated from 

^(R,k) = ^(R,k) + (31) 

poo poo 

/ G+(R,r',k)\r',a> dr'X a (k) / < r> , a\if>+(r', a)dr' . 
Jo Jo 

Using 

k) - ^(R, k) = k)(S a (k) - S° a (k)) (32) 

poo poo 

= / G^(R,r)\r, a > Sunitarydr\ a (k) / < r, a\ip^ (r, k) dr. 
Jo Jo 

and the final separable potential strengths 

(poo poo \ — 1 

l-A Q (ife) y y <r,a|G+(r,r , ,ife)|r / ,a>drdr / J A a (ife). (33) 

The Lippmann-Schwinger equation for the full problem solutions are generated with the 
reference potential solutions and its Green function plus the separable potential with the 
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strengths a a {k) 



*+(r, k)=2p+(r,k) + 



(34) 



oo 



+ 



Ga(r, x, k)\x, a > dx a a (E(k)) 



< l/,«l^ + (l/,«) dy 




Several options of separable potential are available. In principle any rank is admitted but 
this requires to determine the strengths X a (k) from data at several energies around the mean 
energy E(k). The rank > 1 option can be succesfull for single channel cases but suffers from 
lack of energy dependence in data for coupled channels. As we may use any radial form 
factor we found rank one separable potentials in any case the better option and they are: 
a.) Harmonic oscillator radial wave functions <p a (r,n = 1,Hlu); b.) Normalized Gaussians 
N exp — (r — r ) 2 / a o with r = 0.3 — 0.8 fm, a = 0.1 fm, with the option to put the reference 
potential equal to any value for r < r , simulates a smooth boundary condition around r ; 
c.) A step function 0(r o , a) = 1, <p(r±h, a) = 0.5 and zero else, is very suitable to determine 
the boundary conditions at a chosen radius tq. It serves as surface for the transition from 
the hadronic into the QCD domains. It is the usual boundary condition model. Obviously, 
these examples of form factors are may be replaced by microscopic QCD inspired models 
with a wide range of theoretical inspirations 0]. 

D. Derived Quantities 

The probability density of the full 7i problem 



p(r, a, E) = Trace (r, a, E)^/(r, a, E) 



r 



-2 



(35) 



and the loss of flux, from the continuity equation 



0_ 

Of 



■p(r) + (V-j) =0, 



(36) 



(V- j) := Trace - < ^(a, E)\V(.,r) - V f (r, .)\^(a, E) > . 



(37) 



are interesting to display. 
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III. APPLICATIONS FOR vrN AND vrvr SYSTEMS 



We concentrate in this contribution upon mesonic channels and refer for the NN results 
to the contribution by Andreas Funk. In Fig. 1 we show a sample of the latest SAID 7rN 
continuous phase shift solutions (crosses are single energy solutions) and the optical model 
potentials reproduce exactly the continuous soltions (full lines) 0. The rapid changes with 
energy are signaling resonances which have its correspondences in a strong build-up and 
concentration of probability at short distances ~ 0.25 — 0.6 fm. Similar to the probabibilties 
behaves the loss of flux which is concentrated in the same radial region, left and right 
columns. Best fusion examples are seen in the P33, Di 3 , a highly inelastic resonance is Di 5 , 
and no fusion is seen in P 31 . Some resonances feed strongly into other channels, shown 
with the loss of flux which is highly concenrtrated at the same energy, whereas others are 
predominant elastic channel resonances, seeP33. Similar results are observed in other channel 
and in tttt and KN systems ||. We conclude from our calculations that a significant part 
of resonances, listed in the Data Tables, do not fall into the category of fused systems 
but rather show a strong energy dependence with complex structures in probilities and 
absorption, there appearances are more like molecules. Complex and fused structures, like 
dibaryon resonaces, nir see Fig. 2, and other hybrid excitations, are anticipated for many 
hadronic systems and theoretical and experimental work is in progress||. 
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FIG. 2: A sample of tttt probabilities and loss of currents. 5 { shows the p-meson and 5® the /2 



resonance 
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